Shell structure and pairing for interacting fermions in a trap 
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The shell structures for weakly interacting fermions in harmonic oscillator traps at zero temperature 
undergo several transitions depending on the number of particles in the trap and their interaction 
strength. Calculations of the one and two particle spectra give the pairing gaps, and the many 
particle state is constructed by the seniority scheme. For sufficiently few particles, N^ilO 4 , the 
pairing field exceeds the mean-field splitting of single particle levels at the Fermi surface and the 
greater number of almost degenerate states due to the SU(3) symmetry in the harmonic oscillator 
shell lead to supergaps. Deformation and rotation are also discussed. 

PACS numbers: 03.75.Fi, 21.65.+f, 74.20.Fg, 67.60.-g 



Dilute systems of cold degenerate fermionic atoms have 
been trapped and cooled down to a fraction of the Fermi 
temperature [0. Collective modes are studied and the 
superfluid gap is actively searched for. At low particle 
densities the shell structures in traps are pronounced as 
they are in nuclei and the level degeneracies are impor- 
tant for the size of the gaps which can differ substantially 
from those known from homogeneous systems J2| and sys- 
tems with continuous level densities. The atomic traps 
have the advantage as compared to, e.g., nuclei, that one 
can study almost any number of particles and vary their 
density, interaction strength (by tuning Feshbach reso- 
nances ||), and the number of spin states, e.g. 18 for 40 K 
atoms. The object of this paper is to study the novel shell 
structures of such dilute fermionic systems and to deter- 
mine the number of particles and interaction strength, at 
which the larger SU(3) symmetry and degeneracy of the 
harmonic oscillator potential leads to enhanced gaps. 

We treat a gas of N fermionic atoms of mass m in a 
harmonic oscillator (h.o.) potential at zero temperature 
interacting via a two-body interaction with attractive s- 
wave scattering length a < 0. We shall mainly discuss a 
spherically symmetric trap and a dilute gas (i.e. where 
the density p obeys p|a| 3 <C 1) of particles with two spin 
states with equal population. The Hamiltonian is then 
given by 
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For a large number of particles N at zero temperature 
the Fermi energy is for a non-interacting system 



E F = (n F + 3/2)hw ~ (3AT) 1/3 



(2) 



where np is the h.o. quantum number at the Fermi sur- 
face. The h.o. levels are highly degenerate with states 
having angular momenta I — np, np — 2, 1 or due to 
the Z7(3) symmetry of the 3D spherically symmetric h.o. 
potential. However, interactions split this degeneracy. 
In the Thomas-Fermi (TF) approximation the mean-field 
potential is 
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the Fermi energy Ep = h kp(r)/2m+(l/2)muj 2 r 2 +U(r), 
and the density 

p{r) = 4(r)/37r 2 ~ po (1 - r 2 /R 2 TF f 2 , 



(4) 



inside the cloud r 



< Rtf = a sc\/2np + 3, where 

yjh/rrwj is the oscillator length, and p — 
(2n F ) 3 / 2 /3ir 2 al sc the central density ||. These equa- 
tions constitute the energy density functionals which now 
can be solved self-consistently. The latter expression in 
Eq. (||) is valid for dilute systems where the interactions 
contribute a mean-field weak compared to the confining 
potential. 

The system is dilute when the Fermi energy is large 
compared to the mean- field energy, Ep 3> \U\, or equiv- 



alently kp\a\ 
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o\/a sc <C 1. Dense Fermi 



and Bose systems are studied in Refs. |5(] and |6y. We shall 
be particularly interested in very dilute systems where 
also the h.o. energy exceeds the mean-field potentials, 



fiij \U\, or equivalently 



3/2 



a\/a osc < 1 (see Fig. 3). 



The splitting of the degenerate I — up, rip — 2, orO 
states by the mean-field potential can be calculated per- 
turbatively in the very dilute limit. From the radial h.o. 
wave function 1Z„i (r) for the state with angular momen- 
tum I and n radial nodes in the h.o. shell np = 2n + 1, 
we obtain the single particle energies 



np + ^\ ?iuj = J U{r)\K n i{r)\ 2 r 2 dr 
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where the mean-field and wave function overlap is 
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The latter result is obtained using the WKB approxi- 
mation j?J which compares well to numerical results for 
all I when iif > 1 as shown in Fig. 1. The mean-field 
energies are proportional to the coupling and n 3 J 2 and 
are split like rotational bands with a prefactor that is 
smaller than the average mean-field energy by a factor 
~ 3/16. The relative small splitting reflects the fact that 
the mean-field potential is almost quadratic in r and the 
anharmonic terms therefore small. For attractive inter- 
actions (a < 0) the lowest-lying states have small angular 
momentum, which is opposite to nuclei. 
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FIG. 1. Single particle levels for very dilute systems of 
fermions in h.o. shells uf = 2n + I, The mean-field split- 
ting F(uf,1) vs. tif of multiplets I = riF,nF — 2,..., 1 or 
are calculated numerically from Eq. (^). 



The wave function for a pair of particles each with an- 
gular momentum I adding up to total angular momentum 
L is a sum of time-reversed two-particle states 



b L (n F ,l) = y^{lml - m|£0)V'Zm(ri)V ; z 



(8) 



for M = 0. Here Vim(r) = K n i(r)Y lm (9,<j>) is the 
single particle wave function normalized to unity and 
m = —I,..., I. As the spin wave function is always in 
a singlet state, since otherwise the interactions vanish, 
it has been excluded in Eq. (^|). The corresponding gap 
equal to half the pairing energy between two particles, is 
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for even L and zero otherwise. The gap is largest 
when pairs of particles are coupled to total angular mo- 
ment L = for which the Clebsch-Gordan coefficient is 
(Z0Z0I00) 2 = (2Z + 1)- 1 . 

The overlap integrals can be calculated for ny > 1 
within the WKB approximation M 



A (n F ,Z) 
(\a\/a osc )huj 
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Here, the function C(n F ,l) ~ m(7), to leading logarith- 
mic accuracy. Numerical calculations of the gaps are 
shown in Fig. 2. 




FIG. 2. Pairing energies Ao(n F ,l) in 

units of iT^ 1 ^ 2 huj\a\/a osc calculated from Eq. ([j) (see also 
Eq. (|lo|)). For given uf the levels are / = nF,riF — 2, 1 or 
from top to bottom. 



For two or more particles in a shell, it is known from 
nuclear pairing theory that it is generally energetically 
favorable for the particles to pair in states with L = 
where their spatial overlap is maximum. In the degen- 
erate pairing model with interactions only in the L = 
pair state the pairing energy between n v particles in a 
shell is 

» / ^ ( n v — v)(2fl — n v — v + 2) „ , . 

A(n v ,v) = -^ A An ^A , (11) 

where the seniority v is the number of unpaired parti- 
cles. The number of available states in the multiplct 
is ft — 2(21 + 1) for the multiplet I with pairing gaps 
Ao = Ao(riF, I). In the ground state for even n v all par- 
ticles are paired i.e. v — 0, and the lowest excitations are 
reached by breaking pairs ^| . 

When n F (or I) is sufficiently small, however, the 
mean-field splitting of Eq. (||) is smaller than the pairing 
gaps and the pairing acts between all states in an oscil- 
lator shell and not just in a single / multiplet, i.e. the 
full SU (3) symmetry is effectively restored as compared 
to the SU(2) symmetry of the single I multiplet. We 
shall refer to this enhanced pairing as "S'C/(3)-pairing" or 
"super-pairing" . We shall calculate the gap first by writ- 
ing down the "super" -pair wave function and secondly by 
solving the gap equation. 
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Assuming 5/7(3) symmetry the pairing can be calcu- 
lated variationally with a pair wave function that is a gen- 
eralization of Eq. (0) to a sum over I = np, uf — 2, 1 
or 



J2i ai <Po(nF,l) 



(12) 



The weights ai can be found by a variational method 
with the overlap integrals J TZ 2 l {IZ 2 l , v r 2 dr calculated nu- 
merically. For large rip we find a; ~ \/2l + 1 very 
accurately. The sums then reduce to the h.o. level 
density: £* m |^ m (r)| 2 = (47r)- l £,(2i - 
(l/2)dp(r)/dn,F, and the supergap becomes 
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with the particle density p(r) from Eq. (jij). The gap is 
a factor rip larger than that for a single / value because 
the level degeneracy is larger by rip. 

Turning the mean-field splitting back on, the condi- 
tion for superpairing is that the pairing field 2G ex- 
ceeds the mean- field splitting of the single particle levels, 
Using Eqs. (|), (@) and @ the con- 
dition for superpairing becomes 



G4 



rip , 



(15) 



corresponding to iV^lO 3 trapped particles. For np > 
64/5 we find that only the lower lying states /^,64/5 su- 
perpair. The seniority scheme of Eq. ([ll|) also applies to 
many superpairs replacing Ao = G and the level degen- 
eracy for a full h.o. shell fl — (np + l)(np + 2). 

The transition from pairing within a single (np, I) level 
to many levels can be calculated by solving the gap equa- 
tion 



A(n' F> l') 



21' 
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v 

(Mn'p, /')|4tt— <5 3 (n - r 2 )\Mn F , 0) ■ (16) 
m 

within the same shell, i.e., n' F — np. Here, the single par- 
ticle energies are given by Eq. (^) in the very dilute limit 
and the number of particles is such that the Fermi level 
is positioned at p = t nw .i- When tif^S64/5 the gap equa- 
tion reproduces Eq. ([13]), which is followed by a transi- 
tion region as shown in Fig. 3. When ni?^(64/5) ln(nf), 
equivalent to np^SO, the pairi ng g ap approach that for 
a single I level as given by Eq. jlCj). 

When the interaction becomes sufficiently strong pair- 
ing also takes place between different h.o. shells. The 



pair wave function can still be expressed as a sum over 
time-reversed single particle wave functions |J. For in- 
teraction strengths G\n(np)^?iui Eq. (JsJ) is still a good 
approximation and the generalized gap equation is then 
simply Eq. (|l6|), but also summed over n' F ^ np flTT| |. 
For strong interactions, 2G\n(np )^,7iw, the radial wave 
functions differ substantially from the h.o. ones due 
to the strong pair field ]l0| , |ll[| . The coherence length 
£ ~ RtfHui/ A is smaller than the system size and bulk 
superfluidity sets in (see Fig. 3). The finite system be- 
haves like an irrotational fluid and exhibits collective os- 
cillations with frequencies of order ~ Tiuj 




FIG. 3. Diagram displaying the transitions of the level 
structures and pairing gaps A (in units of Tuj) in traps as a 
function of the number of particles N = n F /3 and the scaled 
interaction strength G/hcj = (32/15n 2 )(2n F ) 1 ^ 2 \a\/a osc . The 
dashed lines separate the dense, dilute and very dilute trapped 
gases. The region between the dash-dotted lines indicates 
the transition between superpairing, A = G, and pairing 
in a single Z-level, A ~ (G/n F ) ln(l). Above the full line 
G = hu/2 ln(2riF) inter-shell pairing enhances the gap A 
above Tilu 



Whereas nuclei with partially filled shells generally arc 
quite deformed, 5 = (r 2 ) — 1 ~ 1/np ||, the atomic 
mean-field does not deform in bulk in a spherically sym- 
metric h.o. trap. The h.o. shell at the Fermi surface 
does not deform either because the energy of the super- 
fluid state is always lower than that of the deformed state, 
and the atoms will be a spherically symmetric fluid. 

Deformed traps such as the oblate ones with u± = 
lui = ll>2^uj3 will only have SU(2) symmetry and cor- 
respondingly lower level degeneracies except when ra- 
tios of the trap frequencies are rational numbers. The 
Fermi energy for the deformed h.o. potential is E F ~ 
npl%ijj±_ + riF^huiz, where np now is the h.o. quantum 
number in two dimensions which has degenerate single 
particle states of angular momentum projected along the 
symmetry axis: m = —np, —up + 2, np. 
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For very deformed traps such that uj±_ <C c^iV -1 / 2 , the 
h.o. potential is effectively two-dimensional, U(r) oc (1 — 
r 2 /Rt F )8(Rtf — r), inside the cloud. Such a quadratic 
mean-field potential does not destroy the SU(2) symme- 
try or split the h.o. shells, except for the states with 
lowest angular momentum whose wave functions extend 
slightly outside Rtf- The pairing gaps can be calculated 
for two particles in shell np Jp| exploiting the SO(2,l) 
symmetry |i~4[ of a two-dimensional trap for particles 
with delta-function interactions 



A 2U (S, M) 



h 2 \a\ 
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1 (S-M)\(S + M)\ 

W [(5-M)!(S±M)!]2 
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where M = -S, -S + 2, S - 2, S, is half the total 
angular momentum projected along the symmetry axis, 
and S = 0, 1, 2, is a quantum number similar to 
seniority. The spatial extent of the particles along the 
symmetry axis enters through a^ 1 = J \ipn F 3 ( r 3)\ 4 dr3 = 
(mw 3 /27r^) 1 / 2 . It is a remarkable feature that the pairing 
energies are independent of np in two dimensions. In the 
supershell np , pair energies of the states with S < np — 1 
equal those of the previous supershell np — 1, recursively. 
Thus, for example, the energy spectrum corresponding to 
Eq. ([It]) always has states with excitation 2fiuj± 1 which 
follows from the SO(2,l) symmetry Q also for strong in- 
teractions. For more than two particles in a shell the low 
lying excitations will be dominated by paired states with 
the seniority quantum number describing the number of 
unpaired particles. 

We now turn to experimental observation of super- 
fluidity and gaps which requires temperatures below 
T<A//c B . For temperatures T ~ 0.2E F as reached in j^] 
this would require strongly interacting atoms (e.g. near 
a Feshbach resonance) with inter-shell pairing so that 
hto^A^Ep = nphui. Alternatively, a cloud of ~ 10 3 
6 Li atoms, which have a very large and negative natural 
scattering length, the supergap is A = G ~ 10~ 2 Ep and 
cooling below that would be required in order to measure 
the gap. 

Superfluid gaps can be determined in experiments with 
cold trapped Fermions from response functions for mul- 
tipole modes [[l2| or from the scissors mode [[ll| in a de- 
formed trap. Alternatively, we suggest rotating the de- 
formed trap slowly around an axis perpendicular to the 
symmetry axis. When A^Shco superfluidity reduces the 
moment of inertia from the rigid-body one to that of an 
irrotational fluid, which is smaller by a factor 8 2 for small 
deformations. As for rapidly rotating nuclei Q several 
interesting phenomena such as superdeformation, back- 
bending, fission, etc. may occur in deformed traps. Op- 
tical lattices in current experiments [jl6|| have few atoms 
in each local trap and we thus expect superpairing which 
favors the insulator vs. conductor state. 

In summary, we have calculated the shell structure for 
a dilute system of interacting fermions as a function of 



the interaction strength and number of particles in a h.o. 
trap. For weak interactions, G^Tilj/2 ln(n^), and few 
particles A^;$10 3 in the trap, the U(3) symmetry of the 
h.o. leads to larger degeneracy and superpairing occurs 
with supergaps A = G. For a larger number of particles 
A^IO 4 — 10 5 , the mean-field energy generally dominates 
and splits the t7(3) symmetry of the h.o. into multiplets 
by an energy that is large compared with the gaps. For 
two dimensional traps, however, the mean-field does not 
split the U (2) degeneracy and the gaps display a remark- 
able regularity as given by Eq. ([l7|). Atoms with more 
than two spins have larger level degeneracies for spin in- 
dependent interactions, and the gaps are correspondingly 
larger as is also the case for uniform density [|]j|] . 

We thank G. Bruun and C. J. Pethick for valuable 
discussions. 
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